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ABSTRACT 

Two  digital  computer  programs  have  been 
formulated  for  predicting  performance  charac¬ 
teristics  of  centrifugal  and  axial  fans  which 
may  be  used  to  provide  lift  air  for  large  sur¬ 
face  effect  vehicles.  Pressure  rise  and  flow 
characteristics  of  forwardly  curved,  radial, 
and  backwardly  curved  centrifugal  fans  can  be 
evaluated.  Axial  fans  having  either  inlet 
flow  prerotators  or  exit  flow  straighteners 
can  also  be  handled  in  the  program.  Sample 
results  are  presented  to  illustrate  the  method 
of  calculation  and  the  nature  of  the  results. 
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NOMENCLATURE 


CENTRIFUGAL  FAN  (EQUATIONS  1  THROUGH  40) 

English 

a  -  Effective  radius  of  interblade  circulation,  ft 
A  -  Area,  ft2 
b  -  Blade  height,  ft 
C^  -  Discharge  coefficient 
d  -  Diameter,  ft 

e  -  Characteristic  roughness  size,  ft 
g  -  Acceleration  of  gravity,  ft/sec3 
gc  -  Gravitational  constant,  lbj^-ft/lb^-sec3 
k  -  Loss  coefficient 
P  -  Pressure,  lb^/ft2 
AP  -  Pressure  loss,  lb^/.Tt2 
Q  -  Volume  flow,  ft3 /sec 
r  -  Radius,  ft 

Rjj  -  Blade  height  ratio,  h1/b2 
Rd  -  Diameter  ratio,  d[/d2 
U  -  Impeller  speed,  ft/sec 
V  -  Velocity,  ft/sec 
w  -  Weight  flow,  ll^/sec 
W  -  Weight,  lb^ 
z  -  Number  of  blades. 
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vii 


Greek 


3  -  Blade  angle  (measured  from  tangent),  degrees 

Ps  -  preswirl  angle,  degrees 
6  -  Inlet  clearance,  ft 

r|  -  Efficiency 
p  -  Density,  lbm/ft3 

Y  -  Pressure  coefficient 

AY  -  Pressure  coefficient 

0  -  Flow  coefficient 

w  -  Rotational  speed,  rad/sec. 

Subscripts  and  Superscripts 

D  -  Diffuser 

e  _  Inlet  eye 

E  -  Modified  Euler  line 

imp  -  Impeller 

in  -  Inlet 

leak  -  Leakage 

m  -  Radial  condition 

max  -  Maximum 

min  -  Minimum 

out  -  Outlet 
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r 

RC 

S 

static 

total 

u 

u' 

V 

0 

1 

2 

3 


Relative  velocity 

Recirculation  loss 

Preswirl  angle 

Static  conditions 

Total  conditions 

Tangential  velocity  component 

Corrected  tangential  velocity  component 

Velocity  pressure 

Impeller  inlet 

Impeller  blade  inlet 

Impeller  blade  outlet 

Downstream  of  impeller  outlet 

Theoretical  Euler  line  equation. 


AXIAL  FAN  (EQUATIONS  4l  THROUGH  95) 

English 

A  -  Area,  fts 

c  _  Blade  chord,  ft 

Cp  -  Drag  coefficient 

CT  -  Lift  coefficient 
u 

D  -  Diameter,  ft 

gc  -  Gravitational  constant,  ll^-f t/lbf -sec3 
j  -  Number  of  blade  sections 
n  _  Fan  speed,  rpm 
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K, 


N 

P 

Q 

r 

s 

V 

X 

cL/cD 

DH/Dtip 


Number  of  blades 
Pressure,  lb^/ft3 
Airflow  rate,  ft3/sec 
Radius,  ft 

Blade  pitch  (spacing),  ft 
Velocity,  ft/sec 
Diameter  ratio,  Dy/Dtip 
Lift-to-drag  ratio 
Hub-to-tip  ratio,  XH  =  DH/Dt^p  . 


Greek 

a  -  Angle  of  attack,  degrees 

-  Blade  angle  of  impeller,  degrees 

y  -  Angle  which  relative  velocity  makes  with  the  plane  of 
rotation,  degrees 

6  -  Angle  which  relative  velocity  makes  with  the  axis  of 

rotation,  degrees 

e  -  Swirl  coefficient 

h  -  Efficiency 

G  -  Camber  angle,  degrees 

p  -  Air  density,  lbjn/ft3 

a  -  Solidity,  c/s 

0  -  Flow  coefficient 

Y  -  pressure  coefficient 

AY  -  Pressure  drop  coefficient 

wr  -  Tangential  velocity  component,  ft/sec  . 
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Subscripts  and  Superscripts 


an  -  Annulus 

ax  -  Axial  direction 

D  -  Diffuser 

des  -  Design  condition 

duct  -  Duct 

H  -  Hub 

i  -  Blade  section 
imp  -  Impeller 

j  -  Total  number  of  blade  sections 

m  -  Midspan  of  blade 

o  _  Zero  lift  condition 

P  -  Prerotator 

PR  -  Profile 

S  -  Straightener 

sec  -  Secondary 

st  -  Stall 

T  -  Total 

th  -  Theoretical 

tip  -  Blade  tip 

*  -  Optimum  ccndition 

1  -  Impeller  inlet  condition 

2  -  Impeller  outlet  condition. 
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INTRODUCTION 


Cushion  air  fan  characteristics  play  an  important  part  in 
the  speed/range  performance,  as  well  as  the  riding  quality  of  a 
surface  effect  vehicle  (SEV).  The  most  desirable  fan  character¬ 
istic  is  one  having  as  flat  a  slope  as  possible,  that  is,  a  small 
change  in  pressure  output  for  a  large  change  in  fan  output  flow 
rate  at  a  given  fan  speed.  The  flatter  the  fan  characteristic, 
the  smoother  the  riding  characteristics  of  the  SEV. 

At  this  time,  only  axial  and  centrifugal  fans  appear  to 
have  future  application  to  large  SEV  craft.1  This  is  due  mainly 
to  their  high  efficiencies  but  also  to  their  low-to-moderate 
weight  and  size,  and  their  proven  feasibility  and  reliability. 
Centrifugal  fans  with  backwardly  curved  blades  have  been  used 
extensively  in  the  relatively  small  SEVs  that  have  been  built  to 
date.  Axial  fans  have  seen  a  very  limited  application  on  present 
SEVs,  mainly  due  to  the  cost  advantage  of  the  simpler  centrifugal 
fan,  which  may  not  be  a  decisive  factor  on  larger  craft. 

In  this  report,  digital  computer  programs  are  outlined  for 
predicting  fan  characteristic  curves  for  centrifugal  and  axial 
fans.  One  objective  of  these  programs  was  to  provide  the  fan 
performance  prediction  in  a  relatively  uncomplicated  manner  so 
as  to  enhance  their  usefulness.  A  reasonably  high  accuracy  can 
be  achieved  with  the  programs  which  mainly  require  only  the  geo¬ 
metric  details  of  the  fan  to  make  the  performance  prediction. 

The  centrifugal  fan  prediction  program  is  based  largely  on  the 
work  of  Osborne3  and  Shipway.3  Additional  modifications  to 
account  for  nonuniform  discharge  velocities  and  recirculation  losses 
were  added.  The  centrifugal  program  will  predict  fan  performance 
for  a  forwardly  curved,  radial,  or  backwardly  curved  centrifugal 
fan.  The  axial  fan  prediction  program  was  the  result  of  work  pre¬ 
sented  by  Wallis4  and  Shipway.5  The  axial  program  predicts  per¬ 
formance  for  an  axial  fan  having  either  exit  flow  s traighteners  or 
inlet  flow  prerotator  using  either  isolated  or  cascade  airfoil  data. 

With  these  programs,  it  is  possible  to  predict  the  perfor¬ 
mance  of  a  given  axial  or  centrifugal  fan  with  a  reasonable  degree 
of  accuracy.  Also,  these  programs  will  be  used  to  study  the 
effects  of  certain  variables  on  the  fan  characteristic  curve  in 
order  to  determine  what  most  significantly  affects  the  fan  char¬ 
acteristic  slope. 


1  Superscripts  refer  to  similarly  numbered  entries  in  the  Technical 
References  at  the  end  of  the  text. 
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A  printout  of  the  present  centrifugal  and  axial  fan  charac¬ 
teristic  prediction  programs  with  sample  data  are  shown  in 
appendixes  A  and  B,  respectively.  Also,  separate  nomenclatures 
are  used  for  the  axial  and  centrifugal  fan,  and  no  attempt  has 
been  made  to  unify  them  since  the  similar  expressions  can  have  a 
totally  different  form  for  the  different  fans.  Appendix  C  pre¬ 
sents  the  details  of  the  centrifugal  fan  weight  estimation  proce¬ 
dure. 

PREDICTION  OF  CENTRIFUGAL  FAN  CHARACTERISTICS 

The  following  method  outlines  a  means  of  predicting  centrif¬ 
ugal  fan  characteristic  curves  for  squirrel-cage-type  fans  which 
are  commonly  used  in  SEVs .  The  prediction  method  follows  similar 
methods  presented  by  Osborne3  and  Shipway3  but  several  modifica¬ 
tions  have  been  made.  These  modifications  include  mainly  a  cor¬ 
rection  for  nonuniform  discharge  velocities  suggested  in  Shepherd6 
and  Wislicenus,7  and  a  recirculation  pressure  loss  proposed  by 
Galvas.8  The  characteristic  pressure-flow  curve  is  derived  in 
both  dimensional  and  nondimensional  form  bj  .onsidering  the  theo¬ 
retical  Euler  line  equation  and  substracting  the  various  losses 
from  it.  The  method  requires  only  the  following  fan  geometry 
and  system  parameters  to  predict  results: 

•  Inlet  and  exit  blade  diameter. 

•  Inlet  and  exit  blade  heights. 

•  Inlet  and  exit  blade  angles. 

•  inlet  eye  diameter. 

•  Number  of  blades. 

•  Exit  flow  area. 

•  Tip  speed. 

A  comparison  of  predicted  and  experimental  characteristics  for  a 
HEBA/B-type  centrifugal  fan  is  made,  and  the  effects  of  varying 
this  fan  geometry  on  the  characteristic  curve  is  studied. 

r  LOW  THROUGH  THE  CENTRIFUGAL  FAN 

By  conservation  of  mass  for  steady  incompressible  flow,  the 
inlet  flow  to  the  fan  will  equal  the  flow  out  of  the  fan.  How¬ 
ever,  the  flow  through  the  impeller  is  usually  greater  than  the 
inlet  flow  because  of  a  leakage  flow  between  the  impeller  and 
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Therefore,  in  the  analysis  of  the  fan  and  its  losses,  the  impel¬ 
ler  should  be  based  on  equation  (2)  while  the  analysis  of  the 
inlet  and  outlet  should  be  based  on  equation  ( 1) .  It  should  be 
noted  that  even  when  Qj_n  and  Qout  equal  zero,  Qimp  can  be  greater 
than  zero.  The  method  of  calculating  will  be  shown  later. 

EULER  LINE 

The  Euler  line  describes  the  maximum  theoretical  pressure 
flow  line  for  the  centrifugal  fan.  It  considers  the  fan  to  have 
an  infinite  number  of  blades  of  infinitesimal  thickness  so  that 
the  entire  airflow  through  the  fan  follows  the  blade  contour 
exactly,  thus  producing  the  maximum  pressure  rise  assuming  no 
sources  of  loss  are  present.  Figure  2  shows  an  outlet  and  inlet 
velocity  diagram  for  a  centrifugal  fan  impeller. 


Item  (a)  -  Outlet  Velocity  Diagram 


Item  (b) 

Inlet  Velocity  Diagram 
with  Preswirl 

FOR  ZERO  PRESWIRL ()8S=0) 
Vm.*V| 


Figure  2 

Centrifugal  Fan  Velocity  Diagrams 


In  item  (a)  of  figure  2,  the  tangential  component  of  the  absolute 
outlet  velocity  Va  is  Vu  ,  while  the  radial  component  is  Vm  and 
the  relative  velocity  issVr  .  The  maximum  work  that  can  be2done 
on  the  air  will  be  equal  tosthe  energy  in  the  air  leaving  the 
impeller,  where,  for  radial  inflow: 

maximum  work  =  Qimp  =  2-  V  Ua  ,  (3) 

l 
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where  W  is  a  weight  flow  rate  of  air.  In  terms  of  pressure  rise 
rather  than  energy,  equation  (3)  may  be  written  as: 


P  =  —  vu  u.  . 

to  g  U2  2 


Using  figure  1  and  item  (a)  of  figure  2,  the  following  expression 


for  Vu  is  obtained: 

3 


V  =  ua  -  vm  cot  32 
2  2 


U2  "  cot  ' 

2  nd„b_  2 

2  2 


(5a) 


where 


^imp 

ndsbs 


(5b) 


Thus,  the  theoretical  Euler  line  equation  becomes 


P_  =  —  U 


P  u2  Qimp 


g r  8  gc  nd2b2 


cot  ft  . 


The  calculations  carried  out  in  terms  of  physical  variables 
may  be  converted  to  dimensionless  form  by  the  use  of  two  princi¬ 
pal  dimensionless  variables: 


pressure  coefficient  ¥  =  Pressure  f 

-2-  Uf 
9c  3 
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flow  coefficient  0  = 


flow 


(8) 


d3  U 

3 


S 


The  pressure  coefficient  characterizes  the  dynamic  condition  of 
the  flow  and  represents  the  pressure  divided  by  the  velocity 
pressure  corresponding  to  the  impeller  tip  velocity.9  The  flow 
coefficient  expresses  the  kinematic  conditions  of  the  flow  and 
represents  the  ratio  of  the  axial  velocity  of  the  flow  based  on 
the  fan  diameter  to  the  tip  velocity  of  the  fan  impeller.9  Hence, 
using  the  above  expressions,  equation  (6)  becomes  (in  nondimen- 
sional  form): 


-  1  -  ^imp  -2*  cot  Pe 


TTb, 


(9) 


where  0imp  is  based  on  Q£mp  of  equation  (2). 

EFFECTS  OF  A  FINITE  NUMBER  OF  BLADES 

The  Euler  line  equation  was  derived  by  assuming  the  air 
exits  the  fan  at  exactly  the  trailing -edge  blade  angle.  However, 
this  is  true  only  in  the  case  of  an  infinite  number  of  blades. 
This  is  particularly  important  for  the  impeller  channel  of  a 
centrifugal  machine  where  even  for  ideal  flow  (that  is,  without 
friction,  turbulence,  and  separation),  effects  combine  to  modify 
the  expression  for  the  theoretical  Euler  line  equation. 

One  cause  of  deviation  of  flow  from  the  theoretical  ideal 
is  interblade  circulation.  Stodola10  has  suggested  that  it  is 
the  inertia  effects  of  the  fluid  particles  between  adjacent 
blades  of  a  radial  flow  impeller  that  must  be  considered.  Due 
to  their  inertia,  some  particles  tend  to  retain  their  orientation 
with  respect  to  fixed  axes,  the  result  being  a  circulatory  motion 
relative  to  the  channel.  This  is  similar  to  a  container  of  fluid 
whirled  with  an  angular  velocity  of  at  the  end  of  an  arm,  as 
shown  in  figure  3,  where  a  particle  of  fluid.  A,  tends  to  remain 
stationary.  However,  relative  to  a  point  on  the  container,  B, 
point  A  appears  to  be  rotating  in  a  direction  opposite  but  equal 
in  magnitude  to  that  of  toe  arm.  At  any  radius  "a"  within  the 
container,  the  particle  has  a  relcitive  velocity  of  -a^. 
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Figure  3 

Interblade  Circulation 


Assuming  that  the  impeller  channel  flow  behaves  similarly  and 
that  this  motion  is  superimposed  on  the  outward  flow  between  the 
blades,  the  effect  is  to  reduce  the  value  of  the  tangential  veloc 
ity  Vu2  to  VuJ ,  where  Vu'  =  Vua  -  auu.  The  effective  radius  of 
the  area  between  the  blades  is  difficult  to  assess,  but  is  gen¬ 
erally  taken  as  half  of  the  perpendicular  distance  between  the 
blade  tangents  at  the  impeller  periphery,9  that  is: 


TTda  sin  g3 
2z 


(10) 


where  z  is  the  number  of  blades. 
Since, 


a  i!. 


TTd?  sin  go 


2US  _  tt  sin  IT 


j 


('•WMMrWWe 


and 


Vu3'  -  vUa 


tt  sin  3. 


U„ 


(12) 


The  net  result  of  interblade  circulation  is  to  reduce  the  theo¬ 
retical  pressure  based  on  the  theoretical  ideal  velocity  diagram. 
However,  this  reduction  of  output  potential  is  the  result  of  the 
behavior  of  an  ideal  fluid  and  does  not  represent  a  loss  in  terms 
of  fluid  and  rotor  energies  because  the  impeller  is  not  called 
upon  to  supply  this  additional  energy.  The  significance  of  the 
interblade  circulation  correction  is  that  by  increasing  the  number 
of  fan  blades  for  a  given  fan,  while  maintaining  constant  tip 
speed,  diameters,  and  blade  angles,  higher  pressure  outputs  are 
possible  since  the  flow  more  closely  follows  the  blade  profile 
angles;  thus,  the  theoretical  Euler  line  equation  is  more  closely 
approached  since  interblade  circulation  effects  are  minimized. 

The  reduction  in  theoretical  Euler  pressure  due  to  inter¬ 
blade  circulation  is  expressed  by: 


APb  =  — £  Ua  am 
9c  8 


TT 

—  sm 
z 


and  nondimensionally : 


(15a) 


«b  -  T  Sln  • 


(13b) 


EFFECTS  OF  NONUNIFORM  DISCHARGE  VELOCITY 

Another  effect  which  causes  the  fan  to  deviate  from  the 
theoretical  Euler  line  is  that  the  radial  velocity  may  not  be 
uniform  along  the  discharge  of  the  impeller.  The  pressure  pro¬ 
duced  with  a  varying  velocity  is  less  than  that  produced  with  a 
uniform  discharge  velocity  giving  the  same  net  flow  rate.  For 
a  centrifugal  fan  having  a  linear  variation  of  velocity  across 
the  depth  of  the  channel,  as  shown  in  iteiji  (a)  of  figure  4, 
Shepherd6  showed  that: 
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and  nondimensionally: 
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Item  (a) 
Linear  Profile 


Item  (b) 

Sinusoidal  Profile 


Figure  4 

Nonuniform  Discharge  Velocity  Distributions 


For  the  case  of  a  sinusoidal  flow  distribution,  as  shown  in  item 
(b)  of  figure  4,  Wislicenus7  showed  that: 


lpnd  ' 


0  u;  °lmp 

gc 


COt  P? 


(15a 
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and  nondimens iona lly: 


0.  d9 

AYnd  =  s 

Trb2 


[H 


cot  3. 


(15b) 


Equations  (15a)  and  (15b)  also  hold  for  a  sinusoidal  distribution 
in  the  peripheral  direction  around  the  impeller  and  the  two 
velocity  patterns  may  be  superimposed  giving  a  correction7  factor 
of  [(rr3/8)s  -  1].  The  causes  of  the  nonuniform  velocity  distri¬ 
bution  can  be  the  manner  in  which  the  flow  is  turned  from  the 
axial  to  the  radial  direction  in  the  inlet,  as  well  as  by  real 
fluid  effects  such  as  friction  and  separation  on  the  walls  of 
the  disk,  shroud,  and  blades. 


The  net  result  of  nonuniform  discharge  velocity  is  to 
reduce  the  maximum  theoretical  pressure  predicted  by  the  Euler 
line  equation.  However,  like  interblade  circulation,  this  reduc¬ 
tion  in  output  potential  does  not  represent  a  loss  because  the 
impeller  is  not  called  upon  to  supply  this  additional  output. 
Radial  velocity  distributions  at  centrifugal  fan  outlets,  shown 
by  Brotherhood,11  suggest  that  both  linear  and  sinusoidal  veloc¬ 
ity  distributions  combined  are  generally  present. 


Although  interblade  circulation  and  nonuniform  velocity 
distribution  do  not  require  or  expend  any  energy,  they  do  affect 
the  fan  efficiency  since  they  represent  lost  potential  which 
causes  the  theoretical  Euler  line  to  be  modified  downward.  The 
modified  Euler  line  equation  on  which  efficiency  is  based  now 
becomes : 


PE  =  P.  “  APb  -  APnd  ,  (16a) 

and  nond imens iona lly : 


VE  =  -  AYb  -  AYnd 


(16b) 


Thus,  minimizing  interblade  circulation  and  nonuniform  velocity 
distribution  effects  should  improve  efficiency,  since  the  theo¬ 
retical  Euler  line  condition  will  be  more  closely  approached. 
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IMPELLER  LOSS 

A  loss  occurs  in  the  impeller  due  to  separation  of  the  fluid 
stream  during  flow  through  the  blade  passage  and  is  defined  as: 


Apimp  ^imp 


•  vrJa  • 

2  g_  1  2 


(17) 


The  value  of  k^p  is  between  0.2  and  0.3  for  flat  plate  blades 
and  on  the  order  of  0.15  for  aerofoil  section  blades.3'  3 

Zero  preswirl  has  been  assumed  in  the  calculations  of 
since  APimp  is  very  small  at  low  flow  rates  and  will  not  affect 
the  final  result,  while  at  high  flow  rates  preswirl  is  zero.8'  3 
By  considering  the  inlet  and  outlet  velocity  diagrams  of  figure  2 
and  putting. 


V, 


!imp 


m 


and  V, 


!imp 


»  nd,bl  “““  ^  "  nd  b  ' 


(18) 


then  equation  (17)  becomes: 


AP.  =  -imP- 
imp  2 


r  1 

gc  imp  ^djbj  sin  ^ 


ndaba  sin 


(19a) 


and  nondimensionally : 


"imp  “  ^ 


/ 0imp  d?  \  _ 1 _  1  1 

l  nb2  J  [RdRb  sin  sin  J 


(19b) 


where  Pd  =  dl /dz  and  Fb  =  b,/^. 
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INLET  LOSSES 


The  inlet  loss  is  divided  into  two  parts.  First,  air  enters 
the  impeller  inlet  section, usua lly  through  a  reducing  section, 
as  shown  in  figure  1.  The  airflow  must  expand  from  the  velocity 
at  the  inlet  eye,  Ve,  to  the  axial  velocity,  VQ,  at  the  impeller 
inlet  (or  exit  of  the  reducing  section)  and  be  turned  from  axial 
to  radial  flow.  The  loss  associated  with  the  expansion  may  be 
expressed  as: 


AP 


“•in. 


m. 


_P_ 

9c 


(20) 


Although  the  inlet  conditions  should  be  based  on  Q^n,  the  calcu¬ 
lations  may  be  simplified  with  minimal  error  by  using  since 

the  leakage  condition  is  not  known  until  the  outlet  conditions 
are  known.  Thus,  equation  (20)  becomes: 


APin 


i 


®imp 


(21a) 


and  nondimens iona lly: 


A^in, 


0! 

imp 


(21b) 


The  value  of  kj,n^  ranges  from  0.2  to  0.5*2'  3 

A  second  loss  occurs  as  a  result  of  expansion  or  contraction 
or  the  air  as  it  enters  the  blade  passages.  The  loss  also  takes 
into  account  the  effect  of  preswirl.  British  work0  suggests  that 
for  low  flow  rate  where  VQsVr  ,  the  preswirl  angle,  0S,  is  on  the 
order  of  45  degrees.  At  higher  flow  rates,  where  VQ>Vr  ,  the 
preswirl  angle  is  zero.  This  loss  is  expressed  by: 
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LEAKAGE  FLOW 


Internal  flow  leakage  occurs  primarily  between  the  impeller 
inlet  and  the  casing  inlet  as  shown  in  figure  1.  Also,  a  leakage 
may  occur  around  the  drive  shaft  where  it  enters  the  casing,  but 
this  is  far  less  serious  than  the  inlet  leakage  and  is  not  con¬ 
sidered.  The  leakage  volume  that  flows  back  into  the  inlet  from 
the  impeller  exit  area  may  be  treated  as  flow  through  an  orifice, 
that  is,  the  leakage  volume  may  be  written: 

Qleak  *  cd  ndi  2--tatlC  90  .  (26a) 

’  o 


where  is  a  discharge  coefficient  on  the  order  of  0.6,  pstatic 
is  the  static  pressure  increase  across  the  fan  impeller,  and  5 
is  the  clearance  between  impeller  and  casing.  Written  nondimen- 
sionally: 


^leak  ~  cd  s  n  5 
aa 


^2? 


static 


(26b) 


where  Ystatic  =  Pstatic/( o/9c)  U1 •  In  general,  the  static  pres¬ 
sure  difference  (Ps)  is  nearly  equal  to  the  fan  total  pressure 
(Ptotal)*  due  to  the  static  depression  caused  by  the  inlet  veloc¬ 
ity  pressure  and,  thus,  Ptotal  may  be  used  in  place  of  Pstatic 
without  significant  error. 

For  the  case  of  Qout  =  0,  a  static  pressure  difference 
exists  across  the  fan  impeller.  Thus,  a  leakage  flow  exists, 
and  there  is  a  flow  through  the  impeller,  Q^mp  =  Qieak'  even 
though  the  flow  downstream  is  zero. 

RECIRCULATION  LOSS 

The  recirculation  pressure  loss  is  due  to  the  leakage  flow 
reentering  the  intake  and  causing  addi.onal  work  to  be  done  to 
overcome  dissipation,  turning,  and  mixii  7  losses.  Osborne3  and 
Shipway3  do  not  consider  a  recirculation  loss,  although  Galvas8 
does  consider  a  somewhat  similar  loss.  The  recirculation  loss 
appeared  to  be  best  expressed  as  a  function  of  the  leakage  flow 
velocity  where: 
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where  kleak  appears  to  lie  between  0.2  and  0,5  based  on  recalcu¬ 
lating  data  presented  by  Shipway3  used  to  compare  to  experimental 
data.  Substituting  equation  (26a)  in  (27)  yields: 


ApPC  _  kleak  ^d  ^static  •  (28a) 

and  nondimensionally : 

=  kieak  CdS  ^static  ’  (28b) 


OUTLET  LOSSES 

There  is  almost  always  an  increase  of  flow  area  from  the 
blade  passage  to  the  volute  casing  which  may  range  up  to  2  1/2 
times.  Thus,  there  is  a  tendency  for  retardation  of  flow 
velocity  with  resultant  eddy  formation.  This  is  expressed  as 
an  expans ion -type  loss  by: 


Apout  =  kout  \  —  [V  -  V3J2  •  (29) 

d  9C 

The  coefficient  kQut  will  vary  with  deviation  from  design  condi¬ 
tions,  but  at  maximum  efficiency  it  is  probably  on  the  order  of 
O.4.2*  3  Also, 


(50) 


where  V3  is  the  average  velocity  at  some  measuring  station  of 
area,  A3 ,  which  may  be  at  the  volute  casing  exit.  If  V.  is  mea¬ 
sured  close  to  the  fan  outlet,  the  case  where  V3  =  Va '  is  obtained, 
and  hence  APout  =  0.  From  the  outlet  velocity  diagrams,  item  (a) 
of  figure  2  and  figure  3*  the  following  is  obtained: 
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In  many  present  SEV  installations,  the  centrifugal  fan  is 
mounted  directly  in  a  large  plenum  rather  than  a  volute.  No 
attempt  is  made  to  recover  the  tangential  component  of  velocity, 
such  as  by  mounting  exit  turning  vanes.  In  this  case,  the  outlet 
pressure  loss  is  due  entirely  to  the  dissipation  of  the  tangen¬ 
tial  component  of  velocity13  and  is  stated  as: 


or  nondimensiona lly : 


(55b) 


where  Vy  is  obtained  from  equation  (12). 
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EFFICIENCY 


The  efficiency  of  a  fan  is  generally  defined  as  the  ratio 
of  useful  output  power  to  the  required  input  power  to  the  impel¬ 
ler.  The  input  power  to  the  impeller  is  equal  to  the  product  of 
the  impeller  volume  flow  rate  and  the  modified  Euler  line  pres¬ 
sure  of  equations  (17)  or  ( 18) .  Thus, 


input  power  =  Qimp*  pe  ,  (34a) 

or  nond imensionally , 

input  power  =  0imp*  YE  .  (34b) 


The  outlet  power  '.s  equal  to  the  product  of  the  outleu  pressure 
times  the  outlet  volume  flow  rate.  The  outlet  volume  flow  rate 
is  equal  to  the  impeller  volume  flow  rate  minus  the  leakage  flow 
rate,  which  flows  back  toward  the  fan  inlet  and  is  calculated 
from  equations  (26a)  or  (26b).  The  outlet  pressure  may  be  based 
cn  either  the  total  or  static  pressure  conditions.  The  total 
pressure  condition  is  obtained  from: 


Ptotal  PE  "  APimp  "  APint  ~  APin2  “  APRC  ~  APout  '  (35a) 

and  nond imensionally: 

’’total  -  TE  -  "imp  -  "in,  *  "in,  '  "rc  ‘  "out  •  (35b) 


T) e  static  pressure  condition  is  obtained  from: 


Pstatic  Ptotal  “  APv  ' 


(36a) 


or 


Y 

static 


=  Y.  .  .  -  AY 
total  v 


(56b) 
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where  a?v  and  AYV  represent  velocity  pressure  and 


APV  =  ~  —  v'3 


(if  static  and  total  pressure 
are  measured  at  fa.i  exit). 


(37a) 


IP  3 

-  -  —  V3 


(if  static  and  total  pressure  are 
measured  downstream  of  fan) ; 

(37b) 


and  nondimens ionally: 


AYv=  o  ~ 
2  Su„ 


t  \  . 

f) 


( if  static  and  total  pressure 

are  measured  at  fan  exit). 


(58a) 


1  /Va 

2  u, 


(if  static  and  total  pressure  are 
measured  downstream  of  fan) . 

(58b) 


The  fan  total  efficiency  may  now  be  expressed  as: 


total 


ptotal  Qjmp  "  Qleak  _  ^total  ^imp  ~  ^leak 


(59a) 


The  static  efficiency  is  expressed  by: 


^static  ^imp  ~  ^leak  ^static  ^imp  “  ^leak 


static 


^imp 


(59b) 
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The  efficiency  of  an  impeller  and  volute  (or  diffuser)  will 
always  be  less  than  that  of  just  the  impeller.  The  difference 
will  be  determined  by  the  efficiency  of  the  volute,  r|D,  where 


hD 


Pg  ~  APout  ~  AYqu'c 

PE  »E 


(40) 


In  some  cases,  fan  characteristics  which  are  presented  by  manu¬ 
facturers  are  based  on  measurements  at  the  exit  of  the  fan 
diffuser  or  volute.  However,  efficiencies  quoted  with  these 
fan  characteristics  often  do  not  include  the  outlet  loss  and 
represent  only  the  higher  fan  impeller  efficiency.  Thus,  when 
comparing  experimental  and  predicted  data  it  is  important  to 
know  how  the  experimental  data  were  obtained. 

COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  RESULTS 

Tne  fan  characteristic  prediction  program  was  used  to  pre¬ 
dict  experimental  performance  for  a  60-inch  outside  diameter  (OD) 
backwardly -curved,  airfoil-bladed,  HEBA/B  centrifugal  fan.  The 
results  are  shown  in  figure  5  along  with  all  the  losses.  Agree¬ 
ment  between  experimental  and  predicted  value  is  quite  good,  but 
this  depends  to  some  extent  on  the  selection  of  loss  coefficients 
which  requires  either  experience  or  experimental  background  to 
make  proper  selections.  For  the  present  case,  the  following  loss 
coefficients  were  used: 


k  •  = 

imp 

0.15 

kini  = 

0.28 

k;  p  •“* 

0.7 

kleak  = 

0.3 

kout  = 

0.4. 
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Also,  the  following  physical  variables  describe  the  HEBA/B  cen¬ 
trifugal  fan  system: 

vtip  =  253  ft/s+ 
dj  =  3.46  ft 

d2  =  5.00  ft 

bj  =  1.62  ft 

b2  =  1.18  ft 

=  25  degrees 

j3a  =  50  degrees 

z  =12  blades 

A3  =  21  ft2  (volute  exit) 

de  =  0.9*  di 

cd  =  0.6 

6  =  1/16  inch  per  foot  of  inside 

diameter  (ID). 


A  sinusoidal  nonuniform  discharge  velocity  is  assumed  to  be 
occurring  at  the  impeller  outlet  since  this  helps  simplify  the 
calculations  and  does  not  appear  to  offer  any  significant  error 
in  results  over  assuming  a  linear  profile  or  combination. 


Of  particular  interest  for  SEVs  is  having  a  fan  with  as 
flat  a  characteristic  as  possible,  that  is,  having  a  large  change 
of  flow  rate  for  a  small  pressure  change  at  a  given  fan  speed. 

To  see  what  effect  the  different  variables  have  on  the  fan  char¬ 
acteristic  curve,  most  of  the  above  variables  for  the  HEBA/B 
centrifugal  fan  were  varied  above  and  below  the  actual  values  in 
the  prediction  program.  Items  (a)  through  (h)  of  figure  6  show 
the  results  where  one  variable  at  a  time  was  varied  with  the 
resultant  static  pressure  and  efficiency  being  plotted  against 
the  outlet  airflow  rate. 


^Abbreviations  used  in  this  text  are  from  the  GPO  Style  Manual, 

1975.  unless  otherwise  noted. 
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Item  (g) 

Effect  of  Number 
of  Blades 


Included  in  the  figures  is  an 
estimate  of  how  the  variable 
change  affects  basic  fan 
impeller  and  hub  weight  as 
compared  to  the  original  HEBA/B 
fan  noted  as  curve  2  in  all 
figures.  The  rational  for  the 
weight  estimates  is  snown  in 
appendix  C.  Item  (aj  of  figure 
6  showed  the  most  significant 
effect  on  character istic  slope. 
When  the  fan  inside  diameter 
approached  the  outside  diameter, 
a  significantly  flatter  char¬ 
acteristic  resulted,  although 
peak  efficiency  was  reduced 
compared  to  the  original  fan. 

The  fan  static  efficiency  and 
pressure  are  both  determined 
at  the  fan  volute  exit.  Vari¬ 
ation  of  the  other  parameters 
in  items  (b)  through  (h)  did 
not  significantly  affect  the 
slope  or  peak  efficiency  but 
could  affect  the  pressure 
flow  output.  In  figure  7. 
both  the  inlet  and  exit  blade 
heights  were  changed  together. 
Increasing  the  fan  blade  height 
allows  greater  airflow  through 
the  fan  at  a  slightly  increased 
pressure  and  a  slightly  more 
favorable  slope  without  signi¬ 
ficantly  affecting  peak  effi¬ 
ciency. 


OUTLET  AIRFLOW  RATE,Q0UT ,FTS/SEC 


Figure  7 

Effects  of  Changing 
Blade  Height 


LIMITATIONS  ON  CENTRIFUGAL  FAN  PREDICTION  PROGRAM 

It  lias  been  found  that  size  affects  efficiency  when  geomet¬ 
rically  similar  fans  of  different  size  are  compared. i3  Increas¬ 
ing  the  fan  diameter  for  a  family  of  similar  fans  will  cause  the 
fan  efficiency  to  increase  due  to  a  decrease  in  relative  rough¬ 
ness  (e/d),  where  e  =  characteristic  roughness  size.  The 
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characteristic  roughness  remains  more  or  less  constant  in  manu¬ 
facturing  processes  so  that  relative  roughness  (e/d)  varies 
inversely  as  the  fan  diameter.  The  degree  of  modification  for 
size  effects  is  small  compared  to  the  accuracy  limits  of  the 
above  method  and  has  been  omitted  to  simplify  the  calculations. 

Increasing  the  number  of  fan  blades  in  the  above  program 
has  the  effect  of  improving  output  flow  and  pressure  as  well  as 
efficiency,  since  interblade  circulation  effects  are  minimized 
to  provide  a  higher  modified  Euler  line.  To  avoid  designs  having 
too  many  or  too  few  blades,  the  pitch -to -chord  ratio  of  the  fan 
is  generally  maintained  between  0.5  and  1.5. 14 

AXIAL  FAN  CHARACTERISTICS  PREDICTION 

The  method  of  predicting  axial  flow  lift  fan  characteristic 
curves  was  based  on  work  by  Wallis4  and  Shipway.6  Using  this 
method,  axial  fans  could  be  designed  with  either  prerotator  or 
straightener  vanes,  as  is  usual  with  surface  effect  vehicles. 

The  program  requires  the  basic  geometry  of  the  fan  and  the 
characteristics  of  the  desired  airfoil  blade  to  be  used.  Both 
isolated  airfoil  and  cascade  data  can  be  used  in  the  program, 
depending  on  the  solidity  value  calculated  for  the  fan. 

The  method  for  predicting  axial  fan  characteristics  is 
essentially  divided  into  three  sections: 

•  Preliminary  design 
s  Detailed  design 

•  Off-design  characteristics. 

in  the  preliminary  design,  the  fan  efficiency  and  other  variables 
based  on  given  flow,  pressure,  and  geometry  are  estimated.  The 
detail  design  will  specify  the  geometric  details,  such  as  inlet 
and  exit  angles  of  the  fan  blade  and  the  prerotator  or  straight¬ 
ener  blades.  The  off-design  characteristics  predict  the  varia¬ 
tion  of  fan  pressure  and  efficiency  with  flow  rate  for  a  given 
fan  speed. 
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PRELIMINARY  DESIGN 


The  preliminary  design  is  necessary  to  calculate  and  esti¬ 
mate  some  of  the  variables  that  will  be  needed  in  the  later 
detailed  design.  The  preliminary  design  estimates  the  total  fan 
efficiency  based  on  the  following  input  design  variables: 

•  Total  pressure,  Pt^s 

•  Flow,  Qdes 

•  Outside  (or  tip)  diameter,  Dt^p 

•  Inside  (or  hub)  diameter,  Djj,  or  hub-to-tip 

ratio,  Xh 

•  Tip  speed,  V^ip/  or  fan  r/min,  n. 

Various  combinations  of  the  above  variables  can  be  tried  at  this 
point  to  determine  the  optimum  fan  design  dictated  normally  by 
the  highest  efficiency.  At  the  end  of  the  preliminary  design, 
the  basic  physical  dimensions  have  been  determined  and  a  more 
refined  design  will  be  performed  to  detail  the  aerodynamic 
design  of  the  fan. 

The  preliminary  design  procedure  requires  the  five  input 
variables  above  and  proceeds  to  calculate  total  fan  efficiency 
for  a  fan  having  one  of  two  options: 

•  Prerotators  (or  inlet  guide  vanes) 

•  Straighteners  (or  exit  guide  vanes). 

If  prerotators  are  used,  the  case  is  denoted  by  having  the 
straightener  swirl  coefficient  equal  to  zero  ( es  =  0) .  If 
straighteners  are  used,  then  the  prerotator  swirl  coefficient 
is  zero  ( eD  =  0).  The  swirl  coefficient  represents  the  ratio  of 
the  tangential  velocity  component  to  the  axial  velocity  component 
at  a  given  radius.  The  design  method  used  does  not  cover  the 
case  of  using  both  prerotators  and  straighteners  at  the  same 
time . 


The  axial  velocity  through  the  fan  for  a  given  flow  is  cal¬ 
culated  by: 
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V 


ax 


Dt ip  , 
n  1 


(41) 


-  %*) 


where  Xjj  =  Djj/Dtip.  The  flow  coefficient  at  the  blade  tip  is 
determined  by: 


0 


tip 


Dtif>  vtip 


(42a) 


at  the  hub, 


0H  =  . 

XH 


(42b) 


and  at  the  midspan  of  the  blade. 


<*tiP 


2  0 


tip 


JH 


+  D 


tip 


1  +  X 


H 


2  D 


tip 


( ^2c) 


The  conditions  at  the  tip,  hub,  and  midspan  of  the  blade,  in 
general,  represent  the  extremes  and  average  conditions  associa¬ 
ted  with  the  fan.  Thus,  conditions  at  these  three  sections  are 
important  and  each  is  required  to  be  known  to  assure  a  proper 
axial  fan  design. 

The  theoretical  total  pressure  coefficient  is  calculated 
in  dimensionless  terms  by: 


th 


vtip 


W) 
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For  the  fan  design,  it  has  been  assumed  that  the  total  pressure 
rise  and  the  axial  velocity  component  remains  constant  along  the 
blade  length  and  that  there  is  no  radial  component  of  flow. 

These  assumptions  fulfill  the  requirement  for  the  free  vortex 
flow  that  the  fan  design  is  based  on. 

The  impeller  swirl  coefficient,  which  is  defined  as  the 
ratio  of  the  tangential-to-axial  velocity  component,  is  obtained 
from:4  /  8 


eimptip 


4  vtipS 


n  Vax3  T>r  (1  -  XH3) 


rth  4ip 


(44a) 


eimpm  eimptip 


2 

( 1  +  xh) 


(44b) 


&imPtip 
eirapH  "  XH 


(44c) 


where  r|T  is  the  total  fan  efficiency  and  is  assumed  to  be  0.8 
for  initial  calculations.  The  swirl  coefficient  is  an  important 
variable  in  that  it  represents  a  measure  of  the  torque  of  the 
impeller,  in  this  case,  and  will  be  used  in  many  subsequent 
equations . 

The  product  of  lift  coefficient  (C^)  and  solidity  ( a  =  C/S) 
at  a  section  on  the  blade  is  given  by:4 


(CLa)  =  2 


G  •  _  _ 

imp 


sin 


v  . 

T  imp 


(*>5) 


where  Y^p  is  the  angle  between  the  resultant  velocity  through 
the  fan  impeller  and  the  plane  of  rotation  and  is  expressed  by:4 
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(if  straighteners 
are  used,  ep  =  0)  , 

(46a) 

cr 

( if  prerotators 
are  used,  =  0)  . 

(46b) 


Figures  8  and  9  show  the  geometric  and  velocity  vector  details 
of  the  blades  which  are  used  for  the  above  and  following  equa¬ 
tions-  If  equation  (45)  is  to  be  solved  at  the  hub,  the  appro¬ 
priate  values  of  0jj,  eimpH,  and  Yimpn  at  t^ie  must  used  in 
the  equations.  Likewise,  if  the  value  at  the  midspan  is  desired, 
the  values  of  0^,  ^imp^  and  Yj_mpm  must  be  used.  In  equations 
to  follow,  the  subscript  'H'  means  values  at  the  hub  and  the 
subscript  'm'  means  values  at  the  midspan. 


AXIS  OF 
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Geometric  Details  of  Blades 
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Figure  9 

Velocity  Vectors  for  Fan  System 
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With  increasing  static  pressure  rise  across  the  fan  impeller, 
there  is  a  sharp  increase  in  drag  just  prior  to  blade  stalling. 
Maximum  efficiency  or  optimum  conditions  should  occur  just  before 
this  rapid  increase  occurs.  The  following  relationships4  express 
the  optimum  lift  coefficient  (CL*): 


The  solidity  at  any  blade  section  may  then  be  obtained  by  divid¬ 
ing  equation  (45)  by  (47).  Thus, 


The  value  of  a  at  the  blade  hub  should  be  checked  to  assure  that 
it  is  within  the  design  limits  of  aH  £  1.5.  If  aH  exceeds  1-5. 
the  preliminary  design  input  variables  should  be  changed  (i.e., 
increasing  and/or  Xh  will  decrease  ajj)  and  the  calculations 

up  to  this  point  repeated. 

The  profile,  secondary,  and  annulus  losses  associated  with 
the  impeller  can  now  be  calculated.  The  values  are  obtained  at 
the  midspan  of  the  blade  since  they  should  represent  the  mean 
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loss  for  the  impeller  as  a  whole, 
profile  drag4/6  is: 


The  loss  of  efficiency  due  to 


A^imppp  cDpr 


'L  sin2  Yimp  TT(1  -  *1?)  ' 


(49) 


where  Yimp  is  obtained  from  equation  (46)  and  Cdpr  is  the  pro¬ 
file  drag  coefficient.  The  loss  of  efficiency  due  to  secondary 
losses  is  obtained  by:4/  6 


AYimpsec  CDsec  0  4 _ 

\h  =  CL  *  sin2  Y imp  *  "(1  -  Xh*)  ' 


(50) 


where  Co  is  the  secondary  drag  coefficient.  Combining  equa¬ 
tions  <*sr  and  (50)  yields : 


'imPPR  +  ATimpsec  CDpR  +  CD« 


sin3  Y imp  n(l  -  XH^) 


(51) 


Shipway6  suggests  that  for  this  case  the  value  of  lift,  -to-drag 
ratio  can  be  set  as  constant  for  preliminary  evaluation  as. 


c - +  c - =  20 -5 

D 

PR  sec 


(52) 


The  magnitude  of  the  annulus  drag  loss  was  assumed  by  Wallis4 
not  to  vary  greatly  from  fan  to  fan  and  hence  he  suggests: 
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imPan 


=  0.02  (for  aerofoil-shaped  blades) 


=  0.03  (for  constant  thickness  plate  blades)  . 


(55) 


The  fan  impeller  efficiency  can  then  be  determined  by  adding 
equations  (31)  and  (53)  or. 


«imppR  °'impsec  tTimpan 
'limp'  -  +  -  + 


(54) 


The  loss  in  efficiency  caused  by  prerotator  or  straightener 
vane  pressure  drop  must  be  determined.  if  straighteners  only 
are  used,4 


es  -  ^imp  and  ep  =  0  , 


(55) 


and  also4 


«s_  =  _ cDs*  0 _  4 

Tth  (2.18  -  1.43  c  )  sin2  y  ■  tt(1  -  XH2) 


(  for  0.5  £  eg  <  1.0), 


(56a) 


s  ^  4 

—  -  - 7 -  ’  - -  (for  e  <  0.5)  , 

3  es  sin3  y_  tt(I-X^)  s 


(56b) 


where 


y  =  cot  — 
s  \  2 


*(t) 


(57) 
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The  straightener  and  prerotator  pressure  loss  are  both  calculated 

at  the  midspan  of  the  blade.  CD  is  the  straightener  drag  coef- 

s 

ficient.  It  is  made  up  of  the  sum  of  the  profile  and  secondary 
drag  coefficients  for  the  straightener  blade.  Wallis4  recommends 
the  following  drag  coefficients  be  used: 


Cn  —  C-n  +  Crj 

us  Us-PR  us-sec 


(58a) 


Crj  —  0 . 0l6  , 

Ds-PR 


(58b) 


's -sec 


=  0.018  <%*2 


(58c) 


The  optimum  lift  coefficient  used  in  equation  (58c)  is  obtained 
from:4 


<V  -  2 


r(  1  -  XH')  -  -4  es0 


\  S  “|  1  .  3  7E 


\  40  /  J 


/  r  a  \ 


If  only  prerotators  are  used  then,4 


sp  =  eimp  and  es  =  0  ' 


(60) 


and  also,4 


Ay  CDp*  0 


Yth  sinS  yd  tt(1  -  X H2) 


(for  0.7  *  e  <  1.5)  , 


(6la) 
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or 


4 


v  f> 

3  sin3 
P 


Y 

P 


4 

"( i  -  v ) 


(for  Ep 


<  0.7), 


(61b) 


where 


Y 


P 


(62) 


The  prerotator  drag  coefficient  is  calculated  exactly  as  in 

equation  (58).  The  optimum  lift  coefficient  of  equation  (58c) 
now  becomes,4 


,  /TT(i-xHa)\3 

1  +  \  T% - j 

n(l  -  Xh2  )  +  4  £p0  \a 

40  J  - 


(63) 


The  clearance  losses  can  be  avoided  with  proper  sealing  of  the 
straightener  or  prerotator  blades  to  the  casing.  Thus,  the 
annulus  drag  component  need  not  be  considered  since  it  will  be 
sma 1 1 . 


The  pressure  loss  associated  with  the  diffusion  process 
can  be  obtained  by  assuming  the  diffusion  efficiency.  Provided 
flow  separation  is  avoided,  diffuser  efficiencies  of  nD  =  0.8 
are  obtainable  and  will  be  assumed  in  the  calculations.  For 
an  annular-type  diffuser,  as  shown  in  figure  10,  in  which  0D  = 
3°,  the  diffusion  pressure  loss  coefficient  ( AYp)  is*4 


«D  -  (1  -  rb) 


/ Aimp  \ 
\  ^duct  ) 


r)„  V 


ax 


2  V. 


tip 


a 


a 


where 


Aimp  =  Impeller  annulus  area 
Aduct=  Diffuser  outlet  or  duct  area. 


(64) 
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Figure  10 
Annular  Diffuser 


If  the  duct  diameter  of  the  fan  unit  remains  constant,  the  loss 
coefficient  becomes,4 


(2  -  XHa 


vax 

2  Vtip3 


(65) 


The  efficiency  loss  associated  with  the  diffusion  is  obtained  by 
dividing  equation  (64)  or  (65)  by  (45). 

The  fan  total  efficiency  may  now  be  determined  from: 


(66) 


If  the  nT  calculated  in  equation  (66)  differs  by  more  than 
from  the  value  of  riT,  assumed  initially  in  equation  (44),  repeat 
the  calculation  substituting  the  new  value  of  t|  in  equation  (44) 
and  repeating  all  calculations  up  to  this  point.  If  the  r|T  cal¬ 
culated  is  within  55^  of  the  t>j,  assumed,  then  the  initial  design 
assumptions  are  reasonable  and  design  calculations  may  continue. 
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DETAILED  DESIGN 


At  this  point,  the  preliminary  fan  design  is  completed  and 
the  detailed  fan  design  will  he  undertaken.  Two  design  methods 
are  available,  the  isolated  airfoil  method  or  the  cascade  design 
method.  For  solidity  values  below  0.7 ,  the  isolated  airfoil  data 
are  used,  while  for  solidity  above  1.0  the  cascade  data  are  used. 
In  between,  either  method  may  be  used.  Wallis4  also  recommends 
that  the  straighteners  or  prerotators  be  designed  by  the  cascade 
method  only.  For  the  design  of  all  blade  elements,  free  vortex 
flow  is  assumed,  that  is,  constant  axial  velocity  along  the 
blades . 

In  designing  the  impeller  blading,  an  airfoil  section  must 
be  chosen  such  as  one  of  the  NACA  profiles.  In  the  detailed 
design,  only  the  lift  coefficient  (C^)  and  the  angle  of  attack 
(a)  for  the  design  point  are  required.  When  the  detailed  design 
is  completed,  the  complete  airfoil  section  characteristic  data 
will  be  needed  to  calculate  off-design  performance. 

PREROTATOR  DESIGN 

If  prerotators  are  used,  their  design  must  be  determined 
before  the  impeller  design  can  be  finalized.  To  define  the 
blade  shape,  it  is  divided  into  an  even  number  of  segments  of 
equal  length.  An  even  number  of  segments  will  also  provide 
specifications  of  the  blade  at  a  midspan  location  which  is 
required  for  the  off -design  characteristics.  Thus, 


where  j  is  an  even  number  of  blade  sections  and  i  denotes  the 
diameter  location  from  the  hub  to  the  tip  of  the  blade.  The 
ratio  of  section-to-tip  diameter  is  expressed  by: 


( for  i  =  0  to  j ) 


(68) 
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Likewise,  the  flow  coefficient  at  each  section  diameter  is 
expressed  by: 


0.  . 

0L  =  (for  i  =  0  to  j)  .  (69) 

xi 


The  prerotator  swirl  coefficient  is  then 


e  =  — m-P— (for  i  =  0  to  j)  ,  (70) 

Xi 


and  the  relative  velocity  of  the  flow  coming  off  the  prerotator 
is  at  an  angle  ( 6p)  with  respect  to  the  axis  of  rotation,  defined 


6Pi  ■  tan_1  (ePi)  (for  i  =  0  to  j)  .  (71) 


The  solidity  at  each  section  for  the  prerotator  is  determined 
by:4 


Up.  -  1  (for  epi  <  0.7)  , 


(T2a) 


and 


2  ep.  sin  Yp. 


(for  0.7  s  £pi  <  1.5)  . 


(72b) 


and  Wallis4  assumes  that  with  prerotators,  =  £  satisfies  most 
cases . 

The  camber  angle  ( see  figure  8)  of  the  prerotator  blade  at 
each  section  is  obtained  from:4 


6Pi  ’  1  -  0'.g/ap.  (*°r  i  -  0  to  j) 


(7?) 
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The  blade  chord  is  equal  to. 


nD,  an< 

JTl 

CPi  =  N 

P 


(74) 


where  Np  is  the  number  of  prerotator  blades.  In  general,  it  is 
recommended  that  the  number  of  prerotatcr  or  straightener  blades 
be  greater  than  the  number  of  impeller  blades.5'16  Strieker  and 
Shank1 5  recommend  that  the  number  of  straightener  or  prerotator 
blades  be  two  more  than  the  number  of  impeller  blades  for  opti¬ 
mized  noise  considerations. 


N 


P 


N. 


N. 

imp 


+  2  . 


(75) 


In  order  to  facilitate  manufacturing  of  the  prerotator 
blades,  the  blade  chord  and/or  the  blade  camber  angle  may  be 
kept  constant.  The  values  of  Cp  and  0p  would  be  set  at  the 
values  calculated  for  the  midspan  and  new  crp.,  5p^,  and  e 
coefficients  must  be  calculated  by  using  equations  (74) ,  (73)* 
and  (71)*  respectively. 

IMPELLER  DESIGN 

The  impeller  design  will  depend  on  whether  prerotators  or 
s traighteners  are  used.  The  impeller  blade  is  broken  into  sec¬ 
tions  along  its  diameter  equal  to  those  used  for  the  prerotator 
or  straightener  section.  Thus,  equations  (67)*  (68),  and  (69) 
apply  for  calculating  D^,  ,  and  0j_,  respectively.  if  no  pre¬ 

rotators  are  used,  then  ep^  equals  zero  for  the  design  point  on 
the  impeller.  However,  if  prerotators  only  are  used,  es^  equals 
zero  at  the  design  point  unless  the  prerotator  blades  were  modi¬ 
fied  to  have  constant  Cp  and/or  6p  to  facilitate  manufacturing. 
In  this  case,  there  may  be  a  swirl  component  coming  off  the 
blade  at  the  design  point  and  thus  es^  /  0.  The  swirl  coeffi¬ 
cient  calculated  for  the  modified  prerotator  must  be  subtracted 
from. 


(  es 


GimPtip 


(76) 
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in  order  to  determine  eS£.  If  the  prerotator  blade  has  not  been 
modified,  es^  equals  zero  at  the  design  point. 

The  product  of  lift  coefficient  (Cl)  and  solidity  (aimp^) 
can  be  expressed  by,4 


aimpi  ~  2  ( €s  +  ep)i  s^-n  Yimpi 


(77) 


where  y^mp  is  shown  in  figure  9  and  determined  by,' 


Yimpi  _  tan 


-1 


40i 


n(  1  “  XHa )  -  2  (es  -  ep) 


i  *i. 


(78) 


At  this  point,  the  lift  coefficient  (Cl)  is  selected  from  the 
characteristic  of  the  airfoil  to  be  used.  The  angle  of  attack 
(a)  at  which  it  occurs  is  also  noted.  Then,  dividing  equation 
(77)  the  value  of  Cj,  selected,  the  solidity  (aimpi)  is  deter¬ 

mined.  The  blade  chord  ic  then  obtained  by. 


cimpi 


nDi  aimp^ 


N. 


imp 


(79) 


If  the  number  of  impeller  blades  (Nimp)  is  not  specific,  it  may 
be  approximated  from  the  hub-to-tip  ratio  (Xy)  given  by,16 


6x 


H 


Ki”>P  -  ( 1  -  XH) 


(80) 


The  blade  angle  ( p^)  is  determined  by. 


^  =  ( Yimpi  +  a) 


(81) 
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STRAIGHTENER  DESIGN 

If  the  design  includes  the  straigh toners  only,  then  £p  =  0 
and  thus , 


£imPtip 

Xi 


(82) 


The  solidity  is  obtained  from:4 


cSi  =  -  ( for  es.  s  0.5) 

5 


or 


(85a) 


r 

i 

■ 


' 


Cfo  •  = 


2  e~,  sin  Ve. 

si  si 


(2.18  -  1.45  es;) 


(for  eSi  >  0.5) 


(85b) 


The  camber  angle  of  the  blade  is  determined  from:4 


where 


6Si  =  tan"1  (es.)  . 


(84) 


(85) 


The  straightener  blade  chord  is 


"D. 

l 


-si 


N„ 


(86) 


where  Ns  is  the  number  of  straightener  blades  as  defined  in 
equation  (75)-  As  before,  the  values  of  cg  and/or  0S  can  be 
kept  constant  at  their  midspan  value  to  facilitate  manufacturing. 
However,  new  rTS^,  6S^,  and  es^  values  must  be  determined  from 
equations  (86),  (84),  and  (85),  respectively. 

OFF -DESIGN  PERFORMANCE  PREDICTION 

The  off-design  performance  details  the  pressure-flow  charac¬ 
teristics  of  the  fan  at  a  given  speed  for  other  than  the  design 
flow  condition.  The  basic  assumption  of  the  off-design  theory  is 
that  conditions  at  the  midspan  of  the  blades  are  representative 
of  the  fan  as  a  whole.  Depending  on  the  solidity  values  calcula¬ 
ted  earlier,  it  will  be  known  whether  isolated  or  cascade  airfoil 
data  should  be  used.  The  data  required  from  the  selected  blade 
characteristic  are  the  angle  of  attack  (a)  and  the  corresponding 
lift  coefficient  (cL),  and  lift-to-drag  ratio  (Cl/Cd).  Several 
angles  of  attack  on  both  sides  of  the  design  point  should  be 
selected  to  describe  fully  the  fan  characteristic  curve. 

The  angle  between  the  mean  relative  velocity  to  the  impeller 
and  the  plane  of  rotation  ( Y-j_mp)  is  defined  by: 


Yimp  “  ( ^m  “  a)  • 


(37) 


The  angle  is  then  used  to  calculate  the  swirl  coefficients,4 


,  .  CL  aimp 

(es  +  eP)  =  i  .  v 

^  2  sin  Y 


imp 


cot  Y  j 


1  + 


(cL/cD) 


(88) 


The  value  of  eg  or  tp  is  obtained  from  the  detailed  design  sec¬ 
tion  for  the  midspan  of  the  straightener  or  prerotator,  depending 
on  which  is  used.  If  straighteners  are  used,  the  value  of  £p  is 
assumed  to  equal  zero  which  means  no  preswirl  is  occurring  prior  to 
the  impeller,  which  greatly  simplifies  the  math  involved  without 
significantly  affecting  the  results.4  However,  if  prerotators 
are  used,  es  can  be  nonzero  due  to  the  swirl  component  coming 
off  the  impeller  blades  at  off-design  flows.  This  can  easily 
be  determined  since  Gp  would  already  be  known.  Thus,  the  flow 
coefficient  can  now  be  calculated  at  midspan  from,4 
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—  rr  ft  *  t 


d  -  _ tan  Yimp _ 

nm  y - = - 

[1  +  0-5*  (ea  -  ep)  tan  Ylmp_ 


(39) 


The  theoretical  pressure  coefficient  will  be. 


2  ( es  +  £p) 


2  r|T 


(90) 


The  fan  system  total  efficiency  can  now7  be  expressed  by , 1 


r)T  =  0.98  - 


sin3  Y- 


55-5 


—  +  cot  y . 

C D/  imP 


+  cot  Y . 


Y  /  Y 

P  /  s 
— -  or  - 


Tth  V  Tth 


(91) 


where 


_P_  = 
Yth 


prerotator  loss  =  —7-3 -  - — ■'  —  - 

y  sin3  Yp  \  ep  +  es 


(92a) 


%  .  v  ,  4  CD 

— -•  »  straightener  loss  =  -  •  — 

Y  ,  sin3  Y^,  '-'L 

th  s 


(92b) 
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The  flow  through  the  fan  is  calculated  by. 


TT 


Q  = 


)tip: 


(1  -Xrf)  Vtip  x,„  0m 


(93) 


The  total  pressure  of  the  fan  is. 


P 


T 


(94) 


The  horsepower  required  to  operate  the  fan  is 


PT*  Q 

55C  *  1T 


(95) 


COMPARISON  OF  PREDICTED  RESULTS 

Several  sample  calculations  were  performed  with  the  axial 
fan  characteristic  program  that  has  just  been  summarized.  The 
base  preliminary  design  data  were: 

Dt j  p  =  4.5  ft 

Dhub  =  2.25  ft 
Vtip  =  **00  ft/sec 
PTdes  =  56.0  lb/ft3 
Qdes  =  125°  ft3 /sec. 
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For  the  above  data,  the  fan  characteristics  were  predicted  for 
the  case  with  prerotators  and  with  straighteners,  as  shown  in 
figure  11.  The  figure  shows  that  with  prerotators,  a  flatter 
characteristic  curve  is  obtained  which  is  desirable  for  SEVf.. 
However,  its  peak  efficiency  is  lower  that  for  the  case  with 
straighteners.  These  data  were  determined  from  the  characteris¬ 
tics  of  a  RAF  6E  blade  section  such  as  shown  in  figure  12. 


FAN  TOTAL 
EFFICIENCY, % 


Figure  11 

Comparison  of  Similar  Fan  Designs  with 
Prerotators  and  with  Straighteners 


LIFT  TO  DRAG  RATIO,  C. /C 


t 

I 


Figure  12 

Characteristics  of  RAF  6E  gection 
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To  determine  the  effect  of  blade  characteristics  on  the 
resultant  pressure-flow  characteristics,  a  NACA  blade  profile 
whose  data  are  shown  in  figure  12  was  used  for  the  same  base 
preliminary  design  data  as  the  RAF  6E  section.  Items  (a)  and 
(b)  of  figure  14  show  the  comparative  results  for  the  two  blade 
profiles  with  prerotators  and  with  straighteners ,  respectively. 
The  difference  between  the  predicted  results  for  the  two  blade 
profiles  is  significant.  The  NACA  65  series  profile  provides 
the  flatter,  more  desirable  pressure-flow  characteristic  for 
SEV  use.  However,  its  peak  efficiency  is  slightly  lower. 

In  the  selection  of  the  most  suitable  airfoil  section  for 
the  impeller,  three  main  properties  of  the  section  are  important 
to  examine. 

•  Cl  -  The  lift  coefficient  determines  both  the  solid¬ 
ity  of  the  impeller  and  the  slope  of  the  pressure-flow  charac¬ 
teristic.  High  values  of  CL  are  advantageous  in  providing  both 
low  solidity  and  low  slope. 

•  Cl/Cl  -  The  lift -to -drag  ratio  tends  to  give  high 
efficiencies  for  high  values  of  the  ratio. 

•  ast  -  c.0  -  The  difference  in  angle  of  attack  between 
stall  and  zero  lift  tend  to  yield  increased  operating  range  and 
reduced  characteristic  slope  at  high  values  of  the  difference. 

The  effects  of  fan  diameter  on  the  fan  characteristics  are 
shown  in  items  (a)  and  (b)  of  figure  15  for  an  axial  fan  with 
prerotators  and  with  straighteners,  respectively,  increasing 
the  tip  diameter  (or  decreasing  the  hub-to-tip  ratio)  will 
decrease  the  slope  of  the  fan  characteristic  and  also  lower  peak 
fan  total  efficiency.  The  figures  also  include  an  estimate  of 
how  the  basic  axial  fan  weight  will  vary  for  the  change  in  vari¬ 
ables  compared  to  the  base  fan  design  of  curve  2.  The  axial  fan 
weight  rational  was  obtained  from  reference  ] . 

LIMITATIONS  OF  THE  AXIAL  FAN  PREDICTION  PROGRAM 

In  the  axial  fan  prediction  program,  the  number  of  impeller 
blndes  does  not  affect  the  fan  performance  other  than  to  increase 
or  decrease  chord  length.  This  is  not  absolutely  true  and  exper¬ 
ience  must  be  relied  on  not  to  choose  too  many  or  too  few  blades. 
This  program  is  intended  to  be  used  in  the  range  of  incompres¬ 
sible  flow  for  an  axial  fan  having  either  straighteners  or  pre¬ 
rotators,  but  not  both. 
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Item  (a) 

Comparison  of  Two  Blade  Profiles  with  Prerotators 


FAN  TOTAL 
EFFICIENCY, % 


Item  (b) 

Comparison  of  Two  Blade  Profiles  with  Straighteners 


FAN  TOTAL 
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Effects  of  Blade 


Profile  on  Axial  Fan  Performance 
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Item  (a) 

Effects  of  Tip  Diameter  and  Prerotators 
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Item  (b) 

Effects  of  Tip  Diameter  and  Straighteners 
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AIRFLOW  RATE,  FT3/SEC 


FAN  TOTAL 
EFFICIENCY, % 


Figure  15 

Effects  of  Tip  Diameter  on  Axial  Fan  Performance 


The  axial  fan  prediction  program  was  not  compared  directly 
with  experimental  data  due  to  the  lack  of  suitable  data  for  com¬ 
parison.  Suitable  experimental  data  require  knowing  what  blade 
profile  section  was  used,  a  fact  which  is  not  included  in  any 
obtainable  experimental  data  to  date.  Shipway6  claims  success 
in  applying  a  very  similar  design  method  to  an  actual  fan  design. 

CONCLUSIONS 

The  fan  characteristic  programs  provide  a  relatively  uncom¬ 
plicated  means  of  predicting  centrifugal  and  axial  fan  performance. 
With  these  programs,  such  things  as  the  effects  of  design  changes 
on  the  performance  and  efficiency  can  be  evaluated.  On  the  basis 
of  using  these  programs,  the  following  general  conclusions  were 
drawn : 

•  Any  modification  to  a  fan  which  results  in  a  flatter 
fan  characteristic  slope  also  tends  to  cause  some  reduction  in 
peak  fan-system  efficiency. 

•  A  significantly  flatter  fan  characteristic  slope 
for  a  centrifugal  fan  with  backwardly  curved  blades  can  be 
obtained  by  allowing  the  fan  inside  diameter  to  approach  the 
outside  fan  diameter,  thus  resulting  in  a  narrower  fan  blade. 

•  An  axial  fan  with  straightener  vanes  is  more  effi¬ 
cient  than  a  similar  fan  with  prerotator  vanes.  However,  an 
axial  fan  with  prerotator  vanes  has  a  flatter  fan  characteristic 
slope  than  a  similar  fan  with  straightener  vanes. 

•  The  selection  of  the  axial  fan  blade  profile  or 
blade  characteristics  can  have  a  significant  effect  on  the 
resulting  fan  characteristic  slope. 

•  Decreasing  axial  fan  hub-to-tip  ratio  or  increasing 
tip  diameter  will  cause  a  much  flatter  fan  characteristic  slope 
for  similar  conditions. 

RECOMMENDAT IONS 

•  The  prediction  programs  should  be  compared  against 
any  suitable  experimental  data  that  can  be  obtained  in  order  to 
evaluate  the  programs  more  completely  and  to  understand  the 
selection  of  suitable  design  variables  for  various  cases. 
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•  The  effects  of  various  parameters  and  combinations 
the 'eo f  on  the  fan  characteristic  slope  require  more  study. 

•  The  effect  of  the  number  of  fan  blades  on  the  per¬ 
formance  of  the  fans  should  be  evaluated  separately  in  the  pro¬ 
grams.  Th  i  may  perhaps  be  done  best  by  expressing  it  as  a 
frictional  <pe  of  loss. 
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APPENDIX  A 

CENTRIFUGAL  FAN  CHARACTERISTIC 
PREDICTION  PROGRAM 
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AXIAL  LIFT  FAN  CHARACTERISTIC 
PREDICTION  PROGRAM 
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APPENDIX  C 

CENTRIFUGAL  FAN  WEIGHT  ESTIMATE 

NOMENCLATURE 

-  Inside  blade  height,  ft 

bQ  -  Outside  blade  height,  ft 
d^  -  Inside  fan  diameter,  ft 
dQ  -  Outside  fan  diameter,  ft 
dH  -  Hub  diameter,  ft 

hH  -  Hub  height,  ft 

tj  -  Base  plate  thickness,  ft 

t  -  Upper  shroud  thickness,  ft 

t  -  Averaqe  blade  thickness,  ft 

aavg 

-  Base  plate  weight,  lb 

vr2  -  Upper  shroud  weight,  lb 

W3  -  Blade  weight,  lb 

-  Hub  weight,  lb 

WCF  -  Centrifugal  fan  weight,  lb 
z  -  Number  of  blades 

0^  -  Inlet  blade  angle,  degrees 

pQ  -  Outlet  blade  angle,  degrees 

o  -  Fan  material  density,  lt^/ft3 

1  -  Base  plate 

2  _  Upper  shroud 

5  -  Fan  blades 

4  -  Hub . 


27-404 


C-l 


The  blade  weight: 


W3  = 


zpt. 


cos 


90°  -^4 


0O 


bQ  +  b± 


(5) 


The  hub  weight: 


W, 


TTd^ 

4 


(*0 


The  basic  centrifugal  fan  weight  is  obtained  by  adding  equations 
( 1)  through  (4).  Thus, 


WCF  =  Wl  +  W2  +  W3  +  W4 


(5) 
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